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The escape rate F of the large-spin model described by the Hamiltonian Ti. = —DS^ — HzSz — H^Sx 
is investigated with the help of the mapping onto a particle moving in a double- well potential U{x). 
The transition-state method yields F in the moderate-damping case as a Boltzmann average of the 
quantum transition probabilities. We have shown that the transition from the classical to quantum 
regimes with lowering temperature is of the first order {dV/dT discontinuous at the transition tem- 
perature To) for hx below the phase boundary line hx — hxc{hz), where hx,z = Hx,z/{2SD), and 
of the second order above this line. In the unbiased case {Hz ~ 0) the result is hxc{0) = 1/4, i.e., 
one fourth of the metastability boundary hxm = 1, at which the barrier disappears. In the strongly 
biased limit 5 = 1 - /i^ < 1, one has hxc ~ {2/3f^'^{VS - V2)5^^^ ?^ 0.2345(5^/^ which is about one 
half of the boundary value hxm = (25/3)^/^ ^ 0.54435^/^ The latter case is relevant for experiments 
on small magnetic particles, where the barrier should be lowered to achieve measurable quantum 
escape rates. 

PACS numbers: 75.45. +j, 75.50.Tt 



I. INTRODUCTION 

The two fundamental mechanisms of the escape of a 
particle from a metastable potential well are quantum 
tunneling through the barrier and the classical thermal 
activation over the barrier. The first mechanism is closely 
related to the tunneling level splitting for a particle in ^ 
double-well potential, which was considered by Huncu 
for the ammonia molecule. j-Qther early studies based 
on the WKB approximationaBu treated the ionization 
of atoms in electric -fields, H cold emission |X)f electrons 
from metal surfaces,El and decay of nuclei.l3 Tunneling 
in spin systems was considered much later: Korenblit 
and ShenderEl calculated the ground-state splitting in the 
high-spin rare-earth compounds with the help of a high- 
order perturbation theory, Chudnovskyfl applied the in- 
stanton technique for the Landau-Lifshitz equation to 
calculate the escape rates. The current broad interest 
to the spin-tunneling problem was initiated, however, 
mainly by the application of the instanton method by 
Enz and SchillingE2l and Chudnovsky and Gunther,Eil|aiLd 
the spin- WKB formalism by van Hemmen and Siit6.ll3 

Studying thermally activated escape of a classical par- 
ticle from the metastable minimum of a potential U{x) 
goes back to Kramers,E3 who solved the Fokker-Planck 
equation describing the diffusion of the particle over the 
barrier. For spin systems, the role of thermal agita- 
tion in overcoming energy Jaarriers (described, e.g., bx 
the Stpiier-Wohlfart modelEj) was stressed by Neel.E£l 
Browntj has derived the Fokker-Planck equation for clas- 
sical spin systems and calculated the escape rate in the 



uniaxial model. 

An extensive reference to the thermal activation and 
tunneling of particles can be found in Ref. ^ to mag- 
netization tunneling in Ref. |l^, and to the thermal ac- 
tivation in classical spin systems in Ref. ^ Spin tun- 
neling was reccpfe jobserved in small magnetic particles 
such as ferritincElO'E^ arid,bapiiiiH|je;rritejH3 and in high- 
spin molecules, Mni2AcE3'E§E3EZll23 (see also Refs. |29| ) 
and Feg.E3 

Considering escape at finite temperatures, the first 
idea is to sum the tunneling and thermoactivation es- 
cape rates as stemming from independent channels: F — 
Tq + Fth. Since the thermoactivation rate follows the 
very steep Arrhenius temperature dependence, Fth = 
Fq exp(— A[//T), the transition between quantum and 
classical regimes occurs at the temperature Tq defined 
by F, = Fth(To). Writing F, = Aexp{-B), ignoring 
prefactors and equating the exponents, one obtains the 
estimation 



Tq^"^ = AU/B, 



(1.1) 



where the superscript at Tq says that the ground-state 
tunneling is considered. For T > T^f "^ one has practi- 
cally F = Fth(T), whereas below the transition F = Fg 
is independent of temperature. The transition between 
the two regimes occures on the temperature interval 
AT - Tlf '^/B. Since B oc S, this is much smaller than 

Tq*'"' in the quasiclassical limit, S ^ 1. The simple sce- 
nario above is the prototype for the so-called first-order 
quantum-classical transition of the escape rate, which is 
accompanied by the discontinuity of dV/dT at Tq. 
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It turns out, however, that for common metastable or 
double-well potentials, such as cubic or quartic parabola, 
another scenario is realized. Below the crossover tem- 
perature Tq the particles cross the barrier at the most 
favorable energy level E(T) which goes down from the 
top of the barrier to the bottom of U{x) with lowering 
temperature. Such a regime is called thermally assisted 
tunneling (TAT). The transition from the classical regime 
to TAT is smooth, with no discontinuity oLdT/dT at Tq, 
and the transition temperature is given byEj 

' = ^o/{27t), Qo = V\U"{xsM)\/m. (1.2) 

Here the superscript in Tq denotes the second-order tran- 
sition, Xsad corresponds to the top of the barrier (the 
saddle point of the potential), and ujq is the so-called in- 
stanton frequency. In fact, the above formula is valid in 
the moderate-damping case, in which the Kramers' re- 
sult for the classical escape rate is independent of the 
friction constant rj and coincides with that of the simple 
transition-state theory (TST). In the strong- and weait^ 
damping cases F p oc l/rj and Fq cx rj, respectively,t2l 
and the formula (F2) should be modified (see, e.g., Ref. 
p^ ). The consideration in the moderate-damping case 
is the most simple, and all the results can be obtained 
from the simple quantum TST formula neglecting dissi- 
pation and giving the escape rate as a Boltzmaitftaverage 
of quantum escape rates at different energies.E3E3 The 
quantum-classical transition of the escape rate includ- 
ing the dissipation in the strong- to moderate-damping 
regimes was descdbed with the help of the Caldeira- 
Leggett formalismO in Refs. ^ ^ The re- 

sults show that in the exponential approximation for the 
escape rate only the second derivative cfT/dT^ is discon- 
tinuous at Tq . More-accurate calculations taking into 
account the prefactoil23'E£l show the smoothening of the 

(2) 

transition in the vicinity of Tq due to quantum effects 
and the thermal distribution, so that all the derivatives 
of F(T) behave continuously. 

The terms first- and second-order quantum-classical 
transitions of the escape_rate used aboizjs are due to 
Larkin and Ovchinnikov.a ChudnovskyEJ stressed the 
analogy with the phase transitions and analyzed the gen- 
eral conditions for both types of quantum-classical transi- 
tions. He has shown that for the second-order transition 
the period of oscillations t{E) in the inverted potential 
~U{x) should monotonically increase with the amplitude 
of oscillations, i.e., with the lowering energy E from the 
top of the barrier. If t{E) is nonmonotonic, the first- 
order transition occurs. Quite recently an effective free 
energy F{E) for quantum-classical transitions of the es- 
cape rate of a spin system was writtenji3 the minimiza- 
tion of which determines the escape rate in the exponen- 
tial approximation: F ~ exp(— i^min/?^)- The latter has 

^ + Fq, just as in the Landau 
3 Here a = corresponds to 



the form F = acj) + b(f> + ci 
model of phase transitions.! 
the quantum-classical transition and & = to the bound- 
ary between first- and second-order transitions. 

In a sense, second-order quantum-classical transitions 
of the escape rate are common, whereas the first-order 



ones are exotic and have to be specially looked for. Nev- 
ertheless, a number of systems and processes showing 
first-order transitions are alreeidy known, e.g., a SQUID 
with two JospflhsoD Jiinctions,L3 false vacuum decay in 
field theories ,li30'Lju3-a,nd dcpinning of a massive string 



from a linear defect All these systems have more de- 
grees of freedom than just a particle, thus the search for 
a physical system equivalent to a particle in a potential 
U{x) leading to the first-order transition of the escape 
rate seems to be quite actual. Qualitatively it is clear how 
U{x) should look: The top of the barrier should be rather 
fiat, whereas the bottom should not. In this for 
the rectangular barrier, tunneling just below the top of 
the barrier is unfavorable, the TAT mechanism is sup- 
pressed, and the thermal activation co mpe tes with the 
ground-state tunneling, leading to Eq. ( |l.l| ). Such p-JiSn 
quirement is satisfied, e.g., for the pinning potentialJi3'E2l 
which consists of periodically spaced narrow pits. Here 
the qualitative results can be easily anticipated, at least 
for particles moving in such a potential. However, the 
exact form of the pinning potential is not known. 

A rather simple and experimentally important sys- 
tem satisfying the above requirement is the uniaxial spin 
model in a field described by the Haniiltonian 



H = -DSl - H,^S^ 



H'xSx- (1.3) 

This Hamiltonian can be mapped0iiii0ii onto a 
particle moving in the potential U{x) which has a 
double- well form in the region of field variables h^^z = 



Hx,z/{2SD), S = S 



1/2 satisfying hi^^ 



0^ < rU 

as the original spin model (^_j) in the classical limit .Ell 
The first-order e scap e-rate transition in the unbiased 
{Hz ~ 0) model (|l.3[) for below some critical value 
was found in Ref. |36|. In Ref. ^ the exact value hxc = 1/4 
was obtained with the help of the particle mapping. One 
can get an idea of why the first-order transition should 
occur at small from the fol low ing simple arguments. 
Since tunneling in the model (1.3) is caused entirely by 
the transverse field H^, it becomes very small for 1. 
In this limit the barrier height AU remains finite, and 
the form of the potential near the bottoms should also 
be preserved. The only possibility for the vanishing tun- 
neling rate is that the barrier becomes very thick, with a 
very flat top (see, e.g., Fig. |l]). The latter is just what is 
needed for the first-order quantum-classical escape rate 
transition. 

The aim of this article is to generalize the approach for 
the biased case Hz ^ and to compute the entire phase 
diagram with the boundary line h^dhz) below which the 
transition is first order. We will use the simple damping- 
independent quantum TST formula as the starting point 
for calculations. This requires justification for our spin 
system. It is known that for the model without the trans- 
verse field the thermoactivation escape rate is propor- 
tional to the damping constant, Fq oc 77, for all values of 
r]x3 The models considering hopping over discretCriei!-. 
els for moderate values of S yield the same result EZI'ES 
Such a situation can be thought of as the weak- and 
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strong-damping regimes at the same time. (This is not 
a contradiction, since according to the Landau-Lifshita. 
equation,Ej but not according to the Gilbert equation,E3 
larger values of rj always lead to a faster relaxation.) 
This situation is different from that of a particle, be- 
cause, in terms of polar angles 6 and ip, in the axially 
symmetric case — the spins cross the barrier not 
through the vicinity of a saddle point but through the 
ridge 9 = 6*, where the energy of the spin has a maxi- 
mum. If the transverse field is applied, the spins flow over 
the vicinity of a saddle point 9 = 9* , ip = 0. This brings 
the system closer to the usual situation with particles, 
and the moderate-damnmg_.mg^^ with the damping- 
independent To appears.BElH'BEa The crossover from 
the strong- to moderate-damping regimes was confirmed 
recently by a numerical solution of the Fokker-Planck 
equation for classical spins in the oblique field in Ref. |6^. 
The boundaries of the moderate-damping regime for the 
spin system depend, in addition to rj, on and and 
they are not yet well established. Accordingly, an accu- 
rate description of spin tunneling with dissipation is an 
open problem. For this reason we restrict ourselves in 
this work to the simple damping-independent quantum 
TST approach. 

The structure of the main part of this article is the 
following. In Sec. || the fundamentals concerning the 
particle mapping, WKB approximation and the quan- 
tum TST are reviewed. In Sec. Ill the quantum-classical 



It can be shown that if the coefficients Cm satisfy the 



transition of the escape rate in the unbiased case is stud- 
ied, and the escape rate is calculated in the whole tem- 
perature range including the prefactor. In Sec. IV the 



results are generalized for the biased case. The possibili- 
ties of the experimental observation of different types of 
the escape-rate transition are discussed in Sec. |V|. 



II. PARTICLE MAPPING AND ESCAPE RATE 



A. Particle mapping 



The spji 
mappedL^I-pH^-^ 



em with the Hamiltonian (1.3) can be 
^ onto a particle problem 



Schrodinger equation for t he spin problem (1.3), the coor- 
dinate wave function (2.2) satisfies the Schrodinger equa- 
tion for the particle problem ( p.l[ ) , in the stationary case 
with the same energy leve ls _E„, n = 0, 1, ... , 25. The 
potential U{x) in Eq. (2.1) is given by 



U{x) = S^Du{x), 

u{x) = [hx sinh(a;) — /i^]^ — 2hx cosh(x). (2-5) 

It has a double- well form for and Hz inside the region 
closed by the metastability boundary curvdij 



^2/3 , ^2/3 ^ ^ 



(2.6) 



This simple formula derived 50 years ago has been only 
recently tested in experiment on individual single-domain 
particles.E3 

Finding the extrema of U{x) requires the solution 
of the fourth-order algebraic equation for y = exp(a;) 
and it can be better done numerically. In the unbiased 
case, Hz = 0, the top of the barrier is at a; = and 
it corresponds to the saddle point of the classical spin 
Hamiltonian (1.3). The mininum of U{x) is attained at 
/ix cosh(a;min) = 1- One has 

f/sad = -2S^Dhx, C/min = -^'^^(l + hi), (2.7) 

which yields the barrier height 

AU = L/sad - C/„,in = S^Dil - Kf, (2.8) 

as for the original spin problem. The frequency of small 
oscillations near the bottom of U{x) (the attempt fre- 
quency) 



t^o = VC^"(a;mi„)/m = 2SDsJl-hl (2.9) 

coincides with the ferromagnetic resonance frequency in 
the classical limit 5 — > oo. The instanton frequency ujq 
of Eq. (O) reads 



Cuo = 2SD 



hx). 



(2.10) 



n = ^ + u{x), 

2ni 



dx ' 



m = . (2.1) 

2D ^ ' 



The mapping makes a correspondence between the spin 
wave function s = Sm=-S^ml"^)' where \m) are the 

)rdinate wave 

Cm exp(ma;) 



eigenstates of Sz, and the coordinate wave function 
s 



mt^-s ^{S-m)\{S + m)\ 



where 



f(x) = S[hx cosh(x) — hzx] 



and 



S = S +1/2, hx,z = Hx,z/{2SD). 



(2.2) 

(2.3) 
(2.4) 



It becomes much smaller than ujq for /i^; <C 1, which 
signifies the flat top of the barrier in this limit. The 
expansion of the potential near the top of the barrier in 
the unbiased case has the form 

u{x) = -2hx - hxC^ - hx)x^ + ^ (hx - ^ j a;^ 



(2J1) 



where the change of the sign of the coefficient in the x'^ 
term is responsiblCpfpr the first-order escape-rate tran- 
sition at hx = 1/4.e3 The behavior of u{x) for different 
values of the transverse field is rep resented in Fig. |l|. 

The particle wave functions (2^) imaging those of the 
spin system are proportional to polynomials of power 25* 
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FIG. 1. 



Reduced effective potential for the spin system, 
Eq. (2.5), in the unbiased case. The boundaries between the 
spin states and unphysical states are indicated by horizontal 
dotted lines. 

in y = cxp(a::), thus they can have up to 25* zeros and 
they describe 25 + 1 spin states. The corresponding solu- 
tions of the particle's Schrodinger equation for the potep-i 
tial ( |2.5| ) in the unbiased case were studied by Razavy,EZl 
without a reference to the spin problem. He has found 
explicit analytical solutions corresponding to S = 1/2, 1, 
and 3/2, as well as formally to 5 = 0. The particle prob- 
lem (2.1) possesses, however, an infinite number of states 
with the energy eigenvalues above those of the spin prob- 
lem. These additional unphysical states could, in princi- 
ple, mix together with the true spin states and affect the 
results. This does not happen if the unphysical states are 
much less thermally populated than the spin states near 
the top of the barrier, E — J/gad — —2S^Dhx [here 
is defined by Eq. (2.4) without the tilde]. The criterium 
for neglecting the unphysical states can be formulated, if 
one notices that the boundary between the two groups 
of states corresponds, in the quasiclassical limit, to the 
maximal energy of the spin model (|1.3|), 



(2.12) 



Thus, the unphysical states do not affect the results and 
the analogy with the particle is complete if 



sad 



(2.13) 



This criterium can be violated, if the field hx is small 
and temperature is not low enough. However, we are 
mainly interested here in the temperature range about 
the quantum-classical transition temperature Tq ^ SD. 
For such temperatures the complete analogy with the 
particle is justified by the large spin value S. 



B. Escape rate and level splitting 

The simple quantum transition-state theory postulates 
the escape rate in the quasiclassical case and in the tem- 



perature range T <C Af/ in the plausible form 



1 ^ / 

r = ^2^r„cxpi -■ 



r„, = 



^0 

27r 



W^(-B„), 



T 



2sinh[wo/(2r)]' 



(2.14) 



where ui{En) is the frequency of oscillations at the energy 
level En, W{En) are quantum transition probabilities, 
and Zq is the partition function in the well calculated 
for T <^ AC/ over the low-lying oscillatorlike states with 
En = 1/2)0^0 + ^^111111- Since the transition probability 
W{En) usually increases rapidly with energy E, the sum 
in Eq. ( |2.14| ) extends for not too low temperatures over 
many levels and can be replaced by the integral accord- 
ing to . . . +iW^ dEp{E) . . . with the density of states 
p{E) = l/uj{E)m This leads to 



1 



27rZ, 



Ju, 



dEW{E) exp 



E-U^ 
T 



(2.15) 



For particles C/max — 00 , and the usual expressionU'El is 
recovered. The barrier transparency W{E) is determined 
in the WKB approximation by the imaginary-time action 
(see, e.g., Ref. |68| ) 



S{E) = 2V2m I dx ^JU{x) ~ E, (2.16) 

where xi^2{E) are the turning points for the particle os- 
cillat ing i n the inverted potential —U{x). The factor 2 in 
Eq. ( pl6| ) says that the integral is taken over the whole 
period of oscillations, i.e., the particle crosses the barrier 
twice. The WKB approximation fails for energies near 
the top of the barrier, but the result for W{E) can be 
improved and extended for the energies ab q\r a i hs. bar- 
rier, if the barrier is parabolic near the top.E3c3'E3 The 
result can be conveniently written as 



W{E) 



1 



1 + eMS{E)] ' 



(2.17) 



where S{E) goes linearly through zero for E crossing the 
barrier top level and it should be analytically continued 
into the region E > Usad- In the latter case formula 
( ^.17|) describes quantum reflections for a particle going 
over the barrier, with W{E) slightly lower than 1. In the 
classical limit W{E) ^ 9{E — C/gad), where 9{x) is the 
step function, whereas the partition function Zq given 
by Eq . ( ^.14 ) simplifies to T/loq. Then integration in Eq. 
( p.l5[ ) yields the formula 



r = roexp(-AC//T), 



Zn 



1 — cxp — 



C/n 



Usad 



T 



(2.18) 



If the condi tion ( ^.13 ) is satisfied, the second exponen- 
tial in Eq. (2.18) can be neglected, apd the well-known 
moderate-damping result for particleal3 is recovered. In 
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the opposite case the prefactor Fq in Eq. ( 2.18 ) becomes 
proportional to and vanishes in the limit hx 0. In 
fact, in this case dissipation should be taken into account, 
and the well-known result is Fq oc 77 for all values of the 
damping constant 77. Ej The meaning of this result is that 
in the axially symmetric case precession of a spin around 
the anisotropy axis does not bring it closer to the barrier, 
and the role oiuj{E) as the attempt frequency is lost. In 
this situation spin can overcome the barrier only via the 
diffusion in the energy space. 

The level splitting AEn is related to the barrier trans- 
parency r44^(i?„), and in the WKB approximation it is 
given hy^ 



uj{E^) 



■ exp 



S{E,, 



(2.19) 



For the lowest energy levels this WKB result becomes 
invalid. j-A more accurate consideration of Weiss and 
HaeffnerCHI using the fcmctional-integral technique, 
well as that of ShepardEil improving the usual WKBH' 
scheme shows that for the potentials parabolic near the 
bottom the result above should be multiplied by 



TT (2n-H)"+i/^ 
e 2"e"n! 



(2.20) 



This factor approaches 1 as 1 -I- 1/ (24n) with increasing 
n and it is, in fact, very close to 1 for all n: 1.0750 for 
n = 0, 1.0275 for n = 1, 1.01666 for n = 2, 1.01192 for 
n = 3, etc. Equating the expression in Eq. ( 2.2C| ) to 1 
defines an approximation for n\ which is more accurate 
than the leadi ng te rm of the Stirling formula. The cor- 
rection factor ( 2.20| ), however, does not change much in 
front of the exp onent ially small action term in Eq. (2.19). 

The formula (2.19) has the meaning of the level split- 
ting only if the level s in the two wells are degenerate. 
More generally, Eq. ( 2.19| ) yields the transition ampli- 
tude between the levels on both sides of the barrier, and 
the expression for the escape rate should be sensitive to 
the resonance condition for the lev els in the two wells. 
Thereforej^th^ quantity F„ in Eq. (2.14) should be re- 
placed hyUki 



F„ = 



{AE„ 



■E 



2 

nn' 



(2.21) 



of the precession is completed. This is hardly the case 
for magnetic systems, which is known bpJ-h from exper- 
iments and from theoretical est imat ions .Lj Indeed, reso- 
nant tunneling was observed in experiments on high-spin 
magnetic molecules as well as on underdamped Joseph- 
son junctions.oo In this paper we will, however, ignore 
resonance effects in spin tunneling, which have been con- 
sidered in mo re de tail in Ref. 

Formula (2.14) describes the escape from the 
metastable (left) to the stable (right) well. In the un- 
biased or weakly biased cases the rate of exchange be- 
tween the two wells, which is the observable quantity, 
is the sum of the escape rates from both wells into the 
opposite one (see, e.g., Ref. Thus, in the unbiased 
case all the results for F below should be multiplied by 
2, which will be kept in mind but not done explicitly. 
For low temperatures it is sufficient to introduce a small 
bias field to suppress the backflow from the stable to the 
metastable well and to make the consideration neglecting 
this process absolutely correct. The condition for this is 
i = SH,/T>\. 



III. QUANTUM-CLASSICAL TRANSITION IN 
THE UNBIASED CASE 

A. Level splitting and transition probability 

In the unbiased case , = 0, the imaginary-time ac- 
tion S{E) of Eq. ( |2T^ ) can be reduced to elliptic inte- 
grals, e.g., by the substitution y — cosh(a;). The appro- 
priate formula is available in Ref. and the result is 



S[E) ^ AS[{1 -K)f+Y'^I{o? ,k), 



where 



and 



P = (C/,ad - S)/(f/sad - C/min) 



(3.1) 



(3.2) 



(3.3) 



[see Eq. ( ^.7] )] is the dimensionless energy variable taking 
the values at the top of the barrier and 1 at the bottom 
of the potential. The quantity I{a^ , k) above is given by 



where AS„ is given by Eq. (|IT|), n' are the levels in I = {I + a )K - E + {a - k / a )[n(a ,k) - K] 



the other well and F„„' is the sum of the linewidths of 
the nth and n'th levels. If F„„/ substantially exceeds the 
level spacing u;„/ = E^'+i ~ ^n', the sum in Eq. (2.21) 
can be approximated by the integral and one obtains 



ii-hx)P/f+, k'^f^/u 



(3.4) 



F„ = ^(A^„)V[2c^(S«)]- 



(2.22) 



The lat ter is by virtue of Eq. ( 2.19| ) equivalent to F„ 
in Eq. ( 2.14 ), and the reasoning above can be consid- 
ered as a derivation of it. However, such a situation 
implies that the system is so strongly damped that the 
free precession of the spin decays long before the period 



where K and E are complete elliptic integrals of the 
first and second kind of modulus k and 11(0;^, fc) is the 
complete elliptic integral of the third kind. The latter can 
be expressed through the incomplete elliptic integrals-pf 
the first and second kinds, in our case of < k'^, agI3 



^(a^A:) - K 



a[KE{(3,k) - EF{(3,k)] 



^(l-a2)(fc2-a2) ' 
(3 = arcsin(a/fc). (3.5) 
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The general expression for S{E) simplifies in five lim- 
iting cases. The first case is 1 — /la; <C 1, in which U{x) 
reduces to the quartic parabola [see Eg. (|2.11 ) and Fig. 
Here the parameters a and k of Eq. (3.4)simplify to 



{l-h,)P 



2(1 + 



1 + 



(3.6) 



Expanding Eqs. (3.4) and ( |3.5| ) in powers of ^ 1 one 
obtainsO 

S{E) = (2^/3)[2(l - L)(l + VT^)f'^ 

x[{\ + k'^)E-{\-k^)K]. (3.7) 

For small transverse fields and the energies near the top 
of the barrier, /i^; ~ P <C 1, one can use /+ = P + Ahx 
and /_ = P, which results in fc^ = = P/ {P + Ahx) and 



1 - fc2_ ^ji-jj th&,help of n{P,k) 



E/{1 — fc^) one arrives at the resultli3 

S{E) = 85[P + AK]^''^{K - E 



(3. 



Here for Ahx ^ P <C 1 the modulus k of the ellip tic inte - 
grals is close to 1. In this case using K = ln(4/\/l — fc^) 
and E = well as P = \E\/ [S'^D), one can simphfy 
Eq. (O) tc ' 



S{E) = 4 




(3.9) 



where J7sad is g iven by Eq. (2.7). Now one can rewrite the 
level splitting ( 2.19 ) introducing the unperturbed energy 



levels E„ 



-^,-5 + 1, 



, S and ignoring 



the prefactor in the form 
AS™ ~ {m^/mf^, 



^2S^Kiey8), (3.10) 



which is valid for 1 m <^ S . The latter is the quasi- 
classical limit of the perturbative in result of Refs. |7^, 
^, and ||. 

One can try to continue the perturbative formulas ( |3.9| ) 
and (3.10) to the very top of the barrier, although the 
perturbation theory breaks down there. It can be seen 
that, in contrast to S {E) of Eq. (pM turning to zero at 



E/Usad = 1, Eq. (|3.9| ) turn s to z ero at E/Usad = 



0.92. Accordingly, in Eq. ( |3.10D AP™ ~ 1 at m = nil, 
and not at = 2S^hx, which corresponds to P„i ~ 
—Dm? = C/sad- This can be interpreted as a shift of the 
barrier height by the artifact e^/8 in the perturbative 
formalism. This factor is close to unity and is not very 
important. However, the inaccuracy of the perturbative 
method manifests itself much stronger in the quantum- 
classical transition of the escape rate. The perturbative 
result of Ref. ^ for the boundary between the first- and 
second-order transitions is hxc ~ 0.13, which-is about a 
half of the actual boundary value hxc — l/4.c3 

Another limiting case of Eq. (3T) is the case of a small 
transverse field and energy not too close to the top of 
the barrier, 'ihx <C P ^ 1. This is a purely perturbative 



case, and the level splitting fjoi 
by m is given by the formulaEs 



the energy levels labeled 



AP„ 



2P 



{S+\m\)\ 



[(2H-1)!] {S~-\m\)\ 



2D 



2\rn\ 



(3.11) 



According to Eq. ( 2.19| ) the imaginary-time action S{E) 
can be obtained from the above formula as S{E) = 
— 2 ln(7rAP,„/|a;„i_|_i^m|), where the energy variable P — 
w? /S"^ should be introduced. The WKB approximation 
for S {E) i n th is case, which interpolates between Eqs. 
(p.9[) and (3.19), was calculated in Ref. ^with the help 



of the particle mapping and it can be found there. 

The other two cases in which S{E) simplifies are those 
corresponding to the energy near the top of the barrier 
or the bottom of the well. For the paramctrization it is 
convenient to introduce the dimensionless quantities 



« = ACZ/d-o = ^(1 - hxf'V{2hl/^) 



(3.12) 



and 



^U/loo = ^(1 - hxf'/m + hxf''] (3.13) 



[see Eqs. (^), (^, and (|2.10D]. The quantity v, in 
particular, is a rough estimate of the number of levels 
in the well, based on the assumption that the levels re- 
main equidistant up to the top of the barrier. It thus 
measures quantum effects in the system. One can see 
that for S 3> 1 the system can be made more quantum 
by applying a field close to the metastability boundary, 
1 — hx ^ 1. The condition v ^ 1, which also entails 
w 3> 1, should be, however, satisfied for the applicability 
of the quasiclassical method. The quantity v is related, 
as we will see immediately, to the quantum penetrability 
of the barrier near its top and it also measures quan- 
tum effects. For conventinal potentials, such as cubic or 
quartic parabola, v and v can differ only by a numerical 
factor, and it is sufficient to introduce one of them (see, 
e.g., Ref. This is not the case for our spin model for 
small transverse fields, hx 1- 

Near the top of t he b arrier, P <C 1, the result for S{E) 
following from Eq. (3.1) or directly from the integral ex- 
pression ( 2.16| ) reads 



S{E) = 27rw[P + 6p2 



with 



1 

ih^r 



64 



0(P4)] 



(3.14) 



1 



1 



ihx ^Ghx 



(3.15) 



One can see that b changes sign at hx — 1/4, whereas c > 
for all hx- For small hx the c oeffic ients b, c, etc., become 
large, which means that Eq. ( [3.14|) is only applicable for 
P hx. For hx P I Eq. (|3.8|) can be used. The 
latter gives, however, accurate results only for hx ^ 0.02 
(see Fig. I). 
Near the bottom of the potential one obtains 



6 



s{Eys(u^„) 




FIG. 2. Imaginary-time action S{E) of Eq. for /i^ = 
and different values of and the botto m-ene rgy action 
S{Umin) vs h^. The asymptotes and for 5'(£'), 

as well as S'^'''\Umin) = 43 ln[2 / (eh^)] and Eq. ( po| ) for 
>S(C/min) VS /i are shown by the dotted lines. The result of Eq. 
(p.q) is plotted for hx = 0.01 with the dashed line. 



S{E) - SiU^ 
with 

and 



2{E - U„ 
Wo 



■In 



equoQ 



E-U„ 



In 



1-hl 



\-hl 



(3.16) 



(3.17) 



(3.18) 



q = 8S(1 - hlf'^lhl. 
For /la; <C 1 the bottom-level action S'(C/min) simplifies to 
5(t/„,in) = 45{ln[2/(e/i,)] /i2/4}, (3.19) 

where the first term is the perturbative result, 
5P"'(C/min) (see Refs. |, ^, and Since the correction 
term in Eq. ( |3.19| ) is quadratic in with a small coeffi- 
cient, S'P°''*([/min) is a good approximation to S{Ua\ia) in 
the whole region of the first-order escape-rate transition, 
/la; < 1/4. In the other limiting case one has 



^(f/min) - (45/3)(2e)3/2, 



e = 1 - /la; < 1 (3.20) 



(see, e.g., Ref. |Tl|). The dependences S{E) for different 
values of hx, as well as S{Urain) vs are shown in Fig. 



The instanton period 
dS{E) 



rX2(E) 



-{E)=- 



dx 



dE 



ci(.E) ^U{x) - E 



(3.21) 



can be easier calculated directly than by the differentia- 
tion of Eq. (|3.lD. The result has the form 



x(£)/x((/,ad) 
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FIG. 3. Instanton period t(E) of Eq. ( p.22D for = and 
different values of hx- 



r{E) = [(1 - K)h 



-1/2 



AK/CSD), 



(3.22) 



where bot h /_ 
[see Eqs. (|^ 
P 



and K h ave the same meaning as above 
and (^^)]. Near the top of the barrier, 

0, one has /+ 4/ia;, k ^ Q, K ^ 7r/2, and Eq. 

( |3.22 ) yields the previously known result r = 27r/(I)o with 
luq given by Eq. (2.10). Near the bottom, P ^ 1, one has 
/c ^ 1, and T logarithmically diverges. The dependences 
t{E) for different values of hx are shown in Fig. ^. 

One can see that for hx > 1/4 the period t{E) mono- 
tonically increases with the amplitude of the oscillations 
in the inverted potential —U{x), i.e., with the increase of 
P of Eq. ( |3.3[ ) . In this case the quantum-classical transi- 
tion is second order .Eil For hx < 1/4 the dependence t{E) 
is nonmonotonic, and the transition is first order. Such 
a behavior of t(E) can be easily explained qualitatively. 
For hx > 1/4 the fourth-order term in Eq. (2.11) is pos- 
itive, i.e., U{x) is of the form — + .x*. The inverted 
potential —U{x) is hence of the type x'^ — a;*, which re- 
sults in the increase of r with the oscillation amplitude 
(i.e., with lowering the energy E) and to the second-order 
transition. At hx < 1/4 the anharmonicity of —U{x) has 
the opposite sign, —U{x) ~ x"^ + x'^, which leads to the 
decrease of t when lowering E for energies below the top 
of the barrier. However, with further lowering of E the 
period r begins to increase and diverges logarithmically 
for E approaching the bottom of the potential. 

The period of the real-time oscillations in the potential 
minima can be calculated with the use of the formula dif- 
fering from Eq. ( p. 21 ) by changing sign under the square 
root. The corresponding energy-dependent frequency has 
the form 



uiE) = 2SD[il - hx)f+]"^T:l{2K) 



where 



K = K{r), 



AhxVT^/f. 



(3.23) 



(3.24) 
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FIG. 4. Fr equen cy of the oscillations in the potential wells 
ui[E') of Eq. (3.23) for /iz = and different values of hx- 



and /_!_ is given by Eq. (3.2). Near the bottom, P ^ 1, 
one has /+ — s- l + Zi^, r — > 0, ^ 7''/2, and the frequency 
uj(E) red uces to the previously obtained quantity loq of 
Eq. (2.£). Near the top of the barrier, P ^ 0, one has 
r ^ 1 and uj{E) goes logarthmically to zero. In the case 
Hx = the precession frequency uj{E) can be calculated 
for the original spin model having the energy levels — 
—Dvn? as the energy difference between the neighboring 
levels. The latter is given by w^.m+i — 2Dm, which 
results in lu{E) = 2SDy/P, if we identify P = m'^/S'^. 
The formula (3.23) yields for = the same result with 
S ^ S, which is an immaterial difference for ^ 1. The 
dependences uj{E) for different values of are shown in 
Fig. ^. One can see that they have different types for 
/la; > 1/4 and < 1/4. 

The energy levels £"„ in the wells satisfy the Bohr- 
Sommerfeld quantization rule (see, e.g., Ref. |68|) 

5"'°^'(i;„) EE 2V2^ J dx ^E.a - U{x) = 27r + 

(3.25) 

and can be found numerically. The total number of levels 
in one well is generally given by 

TViev = ^'""^(C/sad)/(27r) (3.26) 

and in the unbiased case has the explicit form 



lev 



25 



arctan 



1 



(3.27) 



In the limit hx ^ Q one has iViev = S. This means that 
one half of the total 25* + 1 = 25 spin levels are in one 
well and the other half are in the other well, thus all the 
levels above the t op of the barrier are unphysical (c f. the 
end of Sec. |IIA| ). In the case 1 - /i^ < 1 Eq. ( |3.27| ) 
simplifies to 



iViev = (45/37r)(l-/i,)3/2. 



(3.28) 



It can be seen that the quantity v of Eq. ( p. 13 ) underes- 



timates the number of levels in the well by factors of 2 
at ft,^ ^ and 16/(37r) « 1.70 at hx 1. 

With the imaginary-time action S{E) and the oscilla- 
tion frequency w{E) having been determined, t he pr ob- 
lem of finding the level splitting AEn of Eq. (2.19) is 



solved. In particular, for energies near the bottom of the 



(|3.16) with 


ting 
form 


i 


^6 



En — Umin = (« + 1/2)^^ 0. T he level split- 
multiplied by the factor ( ^.20 ) takes on the 



AEn^AEoq'^/nl, 



(3.29) 



where q is given by Eq. (3.18) and AEq is the ground- 
state splitting 



AEo 



7ri/2 



exp \ l-hl 



1 



l-hl 



2S 



[l-hlf'hf. (3.30) 



The factor ft,^'^ above signals that the ground-state split- 
ting arises, minimally, in the 25th |Qp4fin of a pertur- 
bation theory in the transverse field.HO'Ej More gener- 
ally, for the excited states one obtains from Eq. W^ ) 
AEn oc h^^ ") JsU^I jj]^ our case 5^1 one can go over 
from 5 = 5+1/2 and hx = hxS/S to 5 and hx. This 
is not an innocent procedure and it could, in principle, 
change the prefactor in Eq. ( |3.30 ). However, the cor- 
responding correction terms cancel each other, and one 
arrives at the same expression without tilde. The lat- 
ter is a particular case of the result obtained by Enz 
and SchillingE2l for a more general biaxial m odel w ith the 
transverse field (see also Ref. 54). Formula (3.30) can be 
rewritten in a more compact form 



5(t/„,in) 



(3.31) 



which looks very similar to the starting formula (2.19). 
The ext ra fa ctor oc S^_^/hx here appears because 

in Eq. ( [3.31 ) the action (3.17) corresponds to the bot- 
tom of the potential, whereas in Eq. ( 2.19 ) the action 
corresponds to the ground-state energy level, Eq. This 
difference due to the zero-point motion between the two 
actions, which is described by Eq. ( 3.16]) is small, but it 
strongly affects the prefactor in Eq. (3.31). 



B. Escape rate in the exponential approximation 

ForTjC AU the dominant contribution to the inte- 
gral ( 2.15| ) is due to the narrow region of energy E where 
the product W{E) exp{—E/T) attains its maximum. Ne- 
glecting the prefactor, one can write 



r^exp(-F„,in/r), 



(3.32) 
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where . 
energy'' 



is the minimal value of the effective "free 



F(PyMJ 



F = E + TS{E) - Un 



(3.33) 



with respect to E. Since with the exponential accuracy 
one can set S{E) = ior E > Usad, the free energy F{E) 
has a downward cusp at E ~ L/gad • Thus the minimum of 
F{E) is attained on the interval [/,nin < E < [/gad- The 
classical regime corresponds, clearly, to the minumum 
just at the top of the barrier, E ~ C/sad- In the quantum 
regime the energy of the minimum shifts down and can 
be found from the condition 



t{E) = l/T, 



(3.34) 



where t{E) = —dS{E)/dE is the period of oscillations 
in the inverted potential — L [(3j) .| Jjhis condition is famil- 
iar from quantum statistics .e3G3'l2I In the instanton lan- 
guage, it determines the trajectory that dominates the 
transition rate at temperature T. 

Above the top of the barrier, P < 0, the effective free 
energy ( p. 33 ) has the trivial form [see Eq. (^^)] 



F{P) = (1 - P)AC/ = + AU, (3.35) 

whereas just below the top of the barrier, < P <C 1, it 
can with the help of Eq. (3.14) be written as 



P(P) = 2TTiT[aP + bP^ + cP^ + 0{P'^)] + AU. (3.36) 
Here a = {T - T^P)/T with 



(3.37) 



see Eqs. (Ilj) and (|2T0| )], the coefficients 6 and c are 
given by Eq. (3.15), and v is given by Eq. (3.12). The. 
analogy with the Landau model of phase transitionsc3 
described by P = acj)^ -\- 60^ -I- c(j)^ + Pq, now becomes 
apparent. The factor a changes sign at the phase transi- 

(2) 

tion temperature T = Tq . The factor b changes sign at 
the field value hx = 1/4 determining the boundary be- 
tween the first- and second-order transitions. The factor 
c remains always positive. The dependence of P on P 
for the entire range of energy is plotted in Fig. [s] with 
the use of the general analytical expression for S{E), Eq. 
(3.1). At hx = 0.3 (Fig. ga) the minimum of P remains 



AU for aU T > T^^\ Below Tp^^^ it continuously shifts 
from the top to the bottom of the potential as temper- 
ature is lowered. This corresponds to the second-order 
transition from thermal activation to thermally assisted 
tunneling, the quantity P playing the role of the order 
parameter. At = 0.1 (Fig. ^), however, there can be 
one or two minima of P, depending on temperature. The 
transition between classical and quantum regimes occurs 
when the two minima have the same free energy, which 
for hx=0.l takes place at Tq = \m^T^^\ 

One can see that the criterion of the second-order 
escape-rate transition is the positiveness of the second 
derivative of the action S{E) and, hence, the effective 
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FIG. 5. Effective "free energy" for the escape rate: (a) - 
hx = 0.3, second-order transition; (b) ~ ~ 0.1, first-order 
transition. 



free energy P(P) defined by Eq. ( 3.33 ). Because of the 
relation 

dr _ d^S _ 1 d^F 
dE ~ '^'dEP ^ ^TdE^ 



(3.38) 



this criterion is equivalent the requirement that the in- 
stantoit«eriod r monotonically increases with decreasing 
energy.H 

The "simple" e stim ation for the crossover tem per ature 
To given by Eq. (1.1) can with the use of Eq. ( p.Sj ) and 
B = S'(C/min) of Eq. (3.17) be explicitly written as 



^(0) _ 



SD 



(l-hx)' 



SD 



\n[2/{ehx)] ' 



hx<^l 



^ ' i-hx«i. 

(0) 



(3.39) 



One can see from Fig. ||b that Tq underestimates the 
crossover temperature. For hx = 0.1 one has Tq"' = 
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FIG. 6. Dependence of the crossover temperature To on the 
transverse field. 




FIG. 7. Dependences of the effective temperature T^s on 
T for different values of the transverse field. 



I.OQIT^^ < Tq. The estimation Tq^^ becomes, however, 
accurate in the limit of small h^. The dependence of the 
crossover temperature To on the transverse field in the 
whole range, < hx < 1, is presented in Fig. ||. The 
temperature dependence of the escape rate can be con- 
veniently written in the form F ~ exp(— AJ7/Teff), where 
the dependence of Tcff = TAU/F^-m on T is presented 
in Fig. for different h^- It can be seen from Fig. 
that the most significant difference between the estima- 
tion Tg"'' and the actual transition temperature Tq arises 
in the limit of a small barrier, that is, for hx I. The 
former is described by the intersection of the dotted Ar- 
rhenius line with the horizontal line corre spon ding to the 
value of rcff(r)/To at T = 0. From Eqs. (|l|) and ( |3.39| ) 



for hx -> 1 one obtains T^ '' /Tf] ' ^ 37r/(8\/2) w 0.833 

The first-order escape-rate transition considered above 
is the transition from thermal activation to thermally as- 
sisted tunneling near the bottom of the potential and not 
directly to the ground-state tunneling. This is due to the 
logarithmic divergence of the instanton period r for the 
energies near J/min- In some field-theoretical models, as, 
e.g., the reduced nonlinear 0(3)-ct model, t approaches 
near the bottom of the potential. Accordingly, the second 
derivative of S{E) and F{E) is negative everywhere, as 
for the rectangular potential for particles. In such a situ- 
ation, as it is clear from Fig. ^b, the minimun of F{E) can 
only be at i!^ = t^ad or E = ?7min- That is, thermal ac- 
tivation competes directly with the ground-state tunnel- 
ing, and the estimation Tq^^ for Tq is exact. Field theories 
showing this extreme case of the first-order escape-rates 
transition were called "type-II theories" , in contrast to 
the "type-I theories" showing a second-order transition. 
In Ref. it was shown that adding a small Skyrme term 
to the reduced nonlinear 0{3)-a model causes r to di- 
verge near the bottom of the potential, with the accord- 
ingly small amplitude. This is, in a sense, similar to the 
situation realized in our spin model for very small hx (see 
Fig. I). 



C. Beyond the exponential approximation 

To obtain the escape rate F wit h prefactor, one should 
perform integration in Eg. (2.15) or, at lower tempera- 
tures, summation in Eq. (2.14). This is of a general in- 
terest in the situation where the second-order quantum- 
classical escape-rate transition is realized in the exponen- 
tial approximation. Here the effect of thermal distribu- 
tion leads to quantum corrections in the classical region 
of temperature and it smoothens the transition to a usual 
crossover without any singularities of r(T). 

In the c ase of a second-order transition, i.e., & > in 
Eq. ( 3.36 ), there are for S ^ 1 four overlapping tem- 
perature ranges where analytical expressions for r(r) 



are available 
[5/(27rw)]i/2 



Above T, 

-1/2 



(2) 



which exactly means a 3> 
l|)], one 



5-1/^ [see Eqs. (|3.36D and 
can neglect the terms with b and c in Eq. (3.14) and ex- 
tend the integpition in Eq. ( 2.15 ) to ±oo. This leads to 
the expressioncJ 



^ ^ ujQ sinh[cJo/(2r)] 
~ 2tt sin[c:;o/(2T)] 

having the asymtote 



' (3-40) 



Wo 

27r 

.(2) 



1 



^0 



A 2 



24T2 



exp 



AU 

't 

AU 



(3.41) 



for T > i.e., for a > 1. For T approaching Tg^^^ 

from above the prefactor in Eq. ( [3.4C| ) diverges because 
of the unlimited contribution from the r ange E < C/sad, 
i.e., P > 0, into the integral in Eq. ( 2.15| ). 

In the temperature range —b <C a <C 1, i.e., across the 
transition region, one can negl ect th e contribution of the 
states with E > C/gad in Eq. (2.15) and write it in the 
form 



F = 



— — / dP cxp[-27:v{aP + bP^ + cP^ 

2ttZo Jq 



(3.42) 



10 



(2) 

where 2ttv = A/7 /Tq . If 6 is jiipt close to zero, one can 
set c = and obtain the resultEj 




exp 



AC/ 



(2) 



46 j 



erfc 




(3.43) 



('2') 

For lower temperatures a< —b, i.e., T < Tq , one can 
no longer use the expansion of S{E) or F{E) near the 
t op o f the barrier, but in this case the integral in Eq. 
(2.15) is dominated by its stationary point. One can 
extend, again, the integration range to ±oo and obtain 
the resultEj 



1 



1 



Zo y/2n\dT/dE\ 



cxp 



F„ 



T 



(3.44) 



and dr/dE are determined by Eqs. (3.33) and 



where F„iin 

( PI and can be calculated numerically. 

It can be see n tha t Eq. (|3.43|) describes the crossover 
between Eqs. ( |3.40D and |3.44D in the narrow region 



Aa ~ [2&/(7rv)] 



around a — 0. In this re- 



gion the erfc function in Eq. (3.43) changes from 2 below 

One can check that the 



(2) (2) 

Tq to small valu es ab ove Tq 



main part of Eq. (3.43), except for the erfc function, is 
the concrete form of Eq. (3.44) in the temperature re- 



(2) I K 

gion just below Tq^ ^ Thus, Eq. (3.44) can be extended 



(2) 

up to Tg by multiplying it by the erfc function. It can- 

(2) 

not however, be ex tended above Tq since in this region 
Tmin in Eq. ( 3.44) doe s not have the same form as its 



equivalent in Eq. (3.43). At the boundary between first- 
and second-order transitions one has 6 = 0, and the term 
with c in Eq. (3.42) should be taken into account. This 
leads to qualitatively similar results; the crossover be- 
tween Eqs. (3.40) and ( ^.44 ) occurs in a narrower region 
Aa ~ c^/^/(27rw)^/'^ ^ S^^^. In the range of transverse 
fields corresponding to the first-order escape-rate transi- 
tion the width of the crossover between the classical and 
quantum regimes is even na rrow er: Aa ^ B^^ ^ S^^ 
[see the discussion after Eq. (1.1)]. 

In the range of the lowest temperatures one should 
take into acco unt quantization of levels in the well and 
use Eq. (2.14) where summation runs near the bottom 
of the potential. Using the osc illato r ener gy le vels £'„ = 
(n -I- 1/2)^0 + C/min and Eqs. ( ^.22] ) and j3.2g| ), one can 
write 



(1 



7r(Ai;o)2 ^ [qe-^o/i2T)^ 



2uJo 



E 

=0 



2n 



(n!)2 



(3.45) 



where the sum i s the modified Bessel function /q. With 
the help of Eq. (3.31) the result can be put into the final 
form 

r - qujoil ~ e-'^''/^)e-^(^""")/o[2ge--«/(2^)]. (3.46) 

Using the asymptotic formula /p (x) = /^/2ttx for 
X 3> 1 one can check that Eq. ( |3.46| ) goes over with 
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FIG. 8. Effective temperature for the escape rate in the 
expon ential app roximation and beyond. Curve 1 uses Eqs. 
( |3.40| ) and ( 3.41: ) and describes the crossover from the classical 
regime with quantum correction s to thermally assisted tun- 
neling (TAT). Curve 2 uses Eq. ( [j.44| ) and describes the TAT 
regi me in the continuous-level approximation. Curve 3 uses 
Eq. (3.45) for the discrete levels and describes the crossover 



from TAT to the ground-state tunneling. 



raising temperature to Eq. ( 3.44 ) with the parameters 
calculated from the action (3.16). The argument of the 
Bessel function in Eq. (3.46) is of order unity for T ^ Tqo, 
where 



Tn, 



Wo 



SDil - hi) 



2\l/2 



2\uq ln[8S'(l - /i2)3/2//i2 



(3.47) 



[cf. Eq. (6.1) of Ref. The temperature Tqo charac- 
terizes the crossover from thermally assisted tunne ling to 
the groun d-sta te tunneling; for T < Tqo Eq. ( |3.4(j ) yield s 
Fo of Eq. ( 2.14 ), multiplied by the correction factor ( 2.2C ) 
squared. For S ^ \ the c rossov er temp eratu re Tqo is 
lower than Tq given by Eq. ( |3^ ) or Eq. because 
of S under the logarithm. 

The results obtained above can be conveniently repre- 
sented in terms of the effective temperature Toff defined 

by 



Wo 

27r 



Q exp 



T 



Wo 

2lT 



exp 



Af/ 



(3.48) 



where the factor Q accounts for the deviation of the pref- 
actor from that of the simple TST. One has, explicitly, 



Teff = TAC//(T„,in-lnQ). 



(3.49) 



The dependence Tcfi(T) is represented in the situation of 
the second-order transition in Fig. ||. One can see that 
all the analytical curves smoothly join each other. The 
escape rate with the accurate prefactor is always higher 
than that in the exponential approximation. For temper- 
atures above the quantum-classical transition this is due 
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FIG. 9. Derivative of the effective temperature for the es- 
cape rate in the exponential approximation and beyond. The 
small circles represent the numerical derivative of the Tgf!{T) 
curve obtained by joining the curves 1, 2, and 3 in the figure 
above. The crossover boun darie s correspond to the argument 
of the erfc function in Eq. (3.43) taking the values ±1. 



to the nonvanishing quantum transparency of the bar- 
rier. At zero temperature the escape rate is higher be- 
cause tunneling occurs from the ground-state level which 
is slightly above the bottom of the potential. The deriva- 
tive dTeff/T is represented in Fig. ^ where the smoothen- 
ing of the second-order escape-rate transition beyond the 
exponential approximation can be clearly seen. 



IV. THE BIASED CASE 

In the general biased case, 0, t he im aginary- 

and real-time actions of Eqs. ( 2.16 ) and ( 3.25| ) can be 
still expr essed in terms of elliptic integrals, since U{x) of 
Eq. (2.5) is proportional to the fourth-order polynom in 
y — exp(a;). We will not do it here because the turning 
points a;i^2(^') for the motion in the potentials ±f7(a;), 
as well as the extrema of U {x) , are given by the solution 
of the fourth-degree algebraic equation and have a cum- 
bersome analytical form. In this section we also restrict 
ourselves to the exponential approximation, since the ef- 
fects associated with the prefactor of the escape rate F do 
not differ qualitatively from those analyzed in the unbi- 
ased case above. We will therefore neglect the difference 
between S and S*, as well as between h^^z and hx,z here. 

It is convenient to start the qualitative analysis from 
the strongly biased case, 5 = I — ^ 1. Here the 



metastability boundary curve of Eq. (2.6) is given by 



In the above expression the term with e^^ which is re- 
sponsible for the formation of the stable minimum of u(a;) 
and is small in the region of interest was dropped. It can 
be checked that the metastable minimum of Eq. ( |4.2| ) dis- 
appears for > h^rn- For = {hxm-K)/hxm < 1 the 
potential u{x) can be approximated by a cubic parabola, 
which will be done in a more general form below. This is 
a standard casejEl in which the second-order escape-rate 
transition takes place. If, on the contrary, the transverse 
field is removed, the barrier height retains a finite value 
but tunneling should disappear, which means that the 
barrier becomes infinitely thick. Indeed, for hx *C hxm 
the third-degree algebr aic equation determining the ex- 
trema of the potential (|4.2|) simplifies to yield 



ymin 
^sad - 



':hx/{2S)[l + hl/{2S% 



(4.3) 



with j/min.sad = exp(x,„in,sad) and 



^sad 



-23/2M'/2[l_/iJ(25/253/2)]. (4.4) 



One can see that in the limit ~* the point Xsnd is 
fixed, Xmin goes to — oo, and the barrier height 



AU ^ S^DS^ 



4\^/^ K 



4 

9 V h 



K 



(4.5) 



remains finite [cf. Eq. (2^)]. As in the unbiased case, 
the very flat top of U (x) favors the first -ord er esc ape- rate 
transition. The solution given by Eqs. (O) and (4^) be- 



comes invalid for hr 



where Xr, 



Xsad- In this 



region the crossover from first- to second-order transition 
is expected. The form of the potential u{x) for different 
values of in the strongly biased case is shown in Fig. 
[l0| . Its remarkable feature is that both Erst- and second- 
order transitions are reahzed for whatever small barrier, 
AU < S^D6^, with the arbitrarily small 6. This is espe- 
cially interesting for the experiments on small magnetic 
particles, S ~ 10^ — 10^, where the barrier should be 
reduced to achieve measurable escape rates. 

In the general biased case simple analytical results for 
the quantum-classical transition temperature Tq can be 
obtained in two limiting cases: for small transverse fields 
and for the fields near the metastability boundary curve. 
In the first case the quantum-classical transition is of 
the first order, and a good estim atio n of the transition 
temperature Tq is given by Eq. ( |1 . l[ ) , where in the ex- 
ponential approximation B — 5'(J7niin)- The analytical 
expression for the bottom-level action S'(L'niin) can be 
found for small and arbitrary from the general for- 
mula ( 2.16| ) . In the most interesting strongly biased case 



(1 - h^J^^I'^ = (2(5/3)^/^ = 0.5443(5^/^. (4.1) this can be done for the arbitrary < h^m with the 



result 



The reduced potential u{x) of Eq. (2.5) simplifies in the 
region of its metastable minimum and maximum to 



i{x) = hl + h4{hx/4:)e-^'' -6e- 



2e^l 



(4.2) 



5(?7„,i„) 9^ 2S^2hxyUn[ - 3V1 - r^in 

+ {rmin + 2)arctanh\/T~^r~7| , 



(4.6) 
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FIG. 10. Reduced potential u(x) of Eq. (|4.2|) for different 
values of /i^ and 5 :s 1 - ft^ < 1. The value /i^ = 0.2345<5^/^ 
corresponds to the boundary between the first- and sec- 
ond-order escape-rate transitions. 



where r^in = ymin/2/min < 1 and y'^^^ = exp{x'^iJ is the 
turning point on the right side of th e barrier correspond- 
ing to the energy Mmin of Eq. (4^). For <C h^m one 
has y'^,^ ^ [^V(2M](1 - hl/S^s^nd r^^^ ^ hi/ 6^ « 1. 
In this case Eq. (pq) simphfies to 



5({/„,i„)-25'^|ln 



27 / it,xm 

2^ 



22 

27 V 



(4.7) 



It is instructive to find 5([/min) for small in the 
whole range of hz with the use of the general perturbatiiis 
expression for the level splitttings in the biased caseE3'E£l 



AE„ 



2D 



[(m' - m- 1)!]2 

I {S + m')\{S - m)\ 
{S - my.{S + m)l 



2D 



(4.8) 



cf. Eq. ( 3.11 )], where m < and m' = — m — Hz/D 



is the matching level on the other side of the barrier. 
For the ground-state level in the metastable well one has 
m = —S and m' — 5(1 — 2hz). Now, in the exponen- 
tial approximation, one can use, according to Eq. (2.19), 

5P"-*(f/min)^^ -2\nAE_s.,S(l-2h^)- 
of the Stirling formula n\ = (n/e) 
result 



Then with the help 
^ one arrives at the 



5P'='-*([/n,in) = 2S{6\n[4Sy{ehx)^] + hzlnhz}. 



(4.9) 

One can check that in the unbiased case {h z — 0, i.e., 
6 = 1) the first (perturbative) term of Eq. (3.19) is re- 



produced, and in the strongly biased case that of Eq. 
(4.7) is recovered. Th e a ccuracy of Eq. (4.£) is not so 
high as that of Eq. (4.7), which contains the impor- 
tant correction term quadratic in hx/h^m- Finally, in 




FIG. 11. The numerically calculated reduced barrier height 
Au = AU / (S^D ) in the strongly biased case. Asymptotes 
(4.5) and (4.1^) shown by the dotted lines reproduce the 
accurate numerical result practically in the whole range of 
hx/h 



the strongly biased case the estimation of the escape- 
rate transition temperature for small transverse fields is 
To ^ T^°^ = AC//SW,ni„) with AU and S{Umin) given 
by Eqs. (U) and ( ^ . 

Near the metastability boundary curve the considera- 
tion begins with the location of the latter from the con- 
ditions u'{xm) — and u"{xm) — 0, which yields the 
equations 



cosh^(a;„) 



1 



sinh^(xm) 



h^ 



where from Eq. ( |2.6| ) follows. Then for the field 

hx — hxni{l ^a;)^ hz — hzmiX ^^)' 



(4.10) 



(4.11) 



where ex,z ^ 1 and {hxm,hzm} are related by Eq. 
one can expand the potential u{x) of Eq. (2.5) near x^n 
in powers of a; — Xm- If hz is not very small, one can 
restrict oneself to the third-order terms and obtain 



u{x) ~ u{Xm) + Ai{x- Xjn) " ^3(2^ " Xm)^ 

with 

A^ - 2eh^/^h^/^ - h'^'^h^'^ 

and 



'^xm'^x 



"'zm '^z- 



(4.12) 
(4.13) 
(4.14) 



The cubic parabola (4.12) is symmetric about Xm and it 
is characterized by the barrier height 

AU = AS^D{2e/?,f^{hx„,hz^f'^ (4.15) 

and the equal real and instanton oscillation frequencies 
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FIG. 12. Imaginary- and real- time actions in the str ongl y 
biased case, Eqs. (1.6) and (1.20). Asymptotes ol Eqs. (1.7), 
(1.21), and (1.22) are shown by the dotted lines. 



uJo=Cjo = 2SD{6ey/'hl/^hli^, 



(4.16) 



where the latter defines the quantum-classical transition 

f 21 A 1 

temperature: To = Tq = a)o/(27r). Equation (4.16) can 
be put into the form 

ujq=Cjq = 257:>(6e)i/4 cot^/^ 61^/(1 + cot^/^ Oh), (4.17) 



sinOji and hzm = cos9h, where, 
(P), hm = (sin2/3 Oh + cos2/3 )-3/2. 



if one sets h^j^ ~ ^r: 
by virtue of Ec 
If, as in Ref. pll, the field changes along the line 9h = 
const, one has = ez = £, which results in e = e. In this 
case Eq. (18) of Ref. |8l| is recovered, where, by defini- 
tion, Wo is twice as small as here. In the strongly biased 
case, however, it is more convenient to have hz fixed and 
swipe hx across the narrow region where the barrier ex- 
ists. Thus one can set hzm. — 1, — (2(5/3)'^/^, ez — 0, 
- ?,2/3^ _ (2(5/3)ea; to obtain 



and e = Km^x 



and 



AC/= (1652^/9)(2e,/3)3/2<s2 



(4.18) 



(4.19) 



Note that at the applicability boundary of the present 
aproach, — 1, one has AC/ — 0.97S^DS'^, which is 
very close to the exact value AU = S^D5^ for ^ h^m 
obtained above (see Fig. ^l]). That is, the barrier height 
scales as AU oc S'^ in the strongly biased case S ^ 1. 
Similarly, the quantum-classical transition temperature 
scales as To cx S. In fact, as will be shown by the numer- 
ical calculations below, this scaling holds not only in the 
strongly biased case, but practically in the whole region 
of hz, excluding that of very small hz- 

As a pendant to Eq. (4.6), one can calculate the real- 
time action given by Eq. ( p.2E| ) for the energy C/gad in the 
strongly-biased case. The result reads 




0.8 h. 10 



FIG. 13. Phase diagram for the quantum-classical es- 
cape-rate transition. The dashed line corresponds to the max - 
imu m of T o as function of h^. The asymptotes of Eqs. (1.24) 
and (1.26) are shown by the dotted lines. 



'""^t/sad) - 2S^/2h^^[ - 3V^ 



(''sad 



2) arctan V?'sad - Ij , 



(4.20) 



where r^ad = J/sad/yLd > 1 and y'^^^ = exp{x'^^^) is the 
turning point on the left side of the meta stable well cor- 
responding to the energy Wgad of Eq. (4^). The top -leve l 
real-time action above determines according to Eq. (3.2(;) 
the number of levels in the well. For h^ ^ h^m one has 
vLd = hx/iU) and r,ad = (2(5)3/2//i^ » 1. Thus in this 
limit 



5^^"'(C^sad) = 2ttSS[1 - 8/{3^/^7T)^hx/hxm] (4.21) 



[cf. Eq. (gj)] and TViev = SS for h.^ = 0. Near the 
metastability boundary the potentia l u{x ) can be ap- 
proximated by the cubic parabola ( 4.121 ). From Eq. 
( P^ one obtains y'^-^ ^ y'^^^ ^ cxp(a;„) = (1 - 



1/3 



1/3 

xrn 



)/h: 

± V^l/(3A3) 



6 /{Shim). More accurately, Xgad,: 



m ± V2eI73 and x^^^ ^ T 

2^2£a;/3. This yields 1 - rmin = 1 - exp(a::min - a^min) - 
and the same for rgad — 1- Expanding now Eqs. 



(4.6) and (4.20) for < 1 results in 



5(t/n,in) = S''^\Usm) = (1655/45) (6e,) 



5/4 



(4.22) 



The beha vior o f S{Urain) and S (Usad) given by Eqs. 
(0) and (^!20|) in the strongly biased case in the whole 
range of hx/h^m is shown in Fig. 

The order of the quantum-classical escape-rate transi- 
tion is determined, as we have seen, by the sign of the 
the coefficient b in the expansion of the imaginary-time 
actio n S{E ) near the top of the barrier [see Eqs. (3.14) 
and ( 3.36| )]. In the general biased case this coefficient is 
given by 
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FIG. 14. 3d plot of To{hx,hz). In the main part of the 
interval < hz < 1 the transition temperature To scales with 

5 = l-h,. 



AU 



(2m)32£:D6 



TT" TJ"" - - (TJ'" 



(4.23) 



[see Eqs. ( |l.2D and (|2j)]. It reduces to that of Eq. ( |3.15| ) 
in the case /iz = where U'J^^ = 0. In the strongly bi- 
ased cas e the potential of the spin system simplifies to 
Eq. (4.2), and there is a simple analytical solution for 
the line hxc(hz) separating the first- and second-order 
escape-rate transitions. To obtain this line, one can cal- 
culate all needed der ivatives of u{x) at the top of the 
barrier from Eq. ( |4.2| ), equate & = in Eq. (4.22), and 
eliminate the terms with h,j. with the help of the relation 
= 2(-2e"^ -f (^e^^) following from u'^^^ = 0. Then 
the condition 6 = yields e^^ = (5/V6, and for = hxc 
one finds 

hxc = (2/3)=^/4(\/3 - \/2)<5^/2 9^ 0.2345,53/2, (4.24) 



I.e., hxc 

^ 0.4308/1 

x7n- Further, for hx 



one has 



u'l^^ = —2(5^(4-^2/3 — 3) resulting in the transition tem- 
perature at the boundary between first-and second-order 
transitions 



3)1/2^, 



(4.25) 



where (4^273 - 3)^/^ = 0.5157 [see Eq. cf. Eq. 

( P7|) ]. 

The maximum of Tq is attained within the range 
of transverse fields corresponding to the second- orde r 
escape-rate transition. In the unbiased case Eq. ( 3.37 ) 

yields evidently rd'"'"'^ = SD /{2Tr) for hx = h^x'"'''^ = 
0.5. In the strongly biased case using the conditions 
du"[hx,x{hx)]/dhx = and u'[hx,x{hx)] = for x{hx) 
corresponding to the maximum of the potential (4.2) one 
obtains e^^ =6/3 and 



for /ii'"""^ = hx„jV2 = 0.384953/2. 



(4.26) 
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FIG. 15. Dependence To(hx) in th e str o ngly biased case (cf. 
Fig. |). Asymptotes based on Eqs. (I.IL (^ )i and (^^, on 
the left side, and on Eqs. (1^) and ( 1.17 ), on the right side, 
are shown by the d ot ted lines. The analytically calculated 
values of Eqs. (4.24)-([l.26|) are shown by diamonds. 



Numerical calculation of the escape rate F in the gen- 
eral biased case in the exponential approximation poses 
no problems. The minimum of the effective free energy 
F{E) of Eq. ( |333| ) with respect to E can be found using 
the imagimary-time action S{E) numerically calculated 
from Eq. (2.16). For any field H one can establish the 
transition temperature Tq, such as for T < Tq the min- 
imum of F(E) no longer corresponds to the top of the 
barrier, as it is in the classical case. Analyzing the depen- 
dence F{E) for different fields allows one to determine the 
order of the quantum-classical escape-rate transition (see 
Fig. H) . The boundary between the first- and the second- 
order transitions can be found the most easily from the 
condition 6 = in Eq. ( 4.23 ), where the derivatives of 
U (x) are calculated at the numerically determined top 
of the barrier. This method yields the same results for 
hxc{hz) as that described above. The resulting phase di- 
agram for the escape-rate transition is shown in Fig. |l^. 
The 3d plot of To{hx, hz) is shown in Fig. |l|. The de- 
pendence To{hx) in the strongly-biased case is given in 
Fig. 111. 



V. DISCUSSION 

We have presented a comprehensive study of the ther- 
mal and quantum decay of a metastable spin state of the 
uniaxial spin system in the arbitrarily directed magnetic 
field. The moderate damping regime has been studied, 
in which the damping does not influence the dynamics 
of the spin system but provides the thermal equilibrium 
with the environment. The method employed is the map- 
ping of a spin system onto the particle in a double-well 
potential, with the subsequent use of the WKB approach. 
The explicit dependence of the escape rate, including the 



15 



prefactor, on temperature, field, and anisotropy constant 
has been worked out and compared with limiting cases 
obtained by others. This calculation shows how formu- 
las describing different regimes join on temperature down 
to lowest temperatures where quantization of levels be- 
comes significant. The crossover from thermally assisted 
tunneling to the ground-state tunneling at T ~ Tqq is 
descr ibed q uantit atively and discussed in detail [see Eqs. 
(3.46) and (3.47)]. The fascinating new feature of this 
analysis is the existence for spin systems of both, first- 
and second-order, transitions from thermal activation to 
thermally assisted tunneling at T = Tq. The kind of 
transition depends on the strength and the direction of 
the magnetic field. We have calculated the boundary in 
the Hx, Hz plane separating the two different regimes. 

The direct analogy with phase transitions exists in 
the limit of a very large spin S. In that limit the de- 
pendence of the transition rate on temperature changes 
abruptly at Tq, pretty much as thermodynamic quanti- 
ties do in the theory of phase transitions. For finite 5, 
both the first- and second-order transitions of the escape 
rate are smeared, similar to the smearing of the phase 
transition in a finite-size system. The reduced width of 
the crossover between the classical and quantum regimes 
AT/Tq is of the order of 5^^/^ for the second-order tran- 
sition and for the first-order transition. 

For a moderately large spin, S ~ 10, one can explore 
in experiment the entire phase diagram shown in Fig. 
p^ . This is the case of molecular magnets like Mni2Ac. 
If greater spins are studied, the external magnetic field 
must be adjusted such that the energy barrier becomes 
small enough to provide a significant tunneling rate. In 
terms of the phase diagram this means that one h as t o 
work close to the metastability boundary line, Eq. (2.6), 
separating the field ranges with and without the barrier. 
According to Fig. |l3|, for a uniaxial spin close to the 
metastability boundary, both, first- and second-order, 
transitions coexist only in the lower right corner of the 
phase diagram. In this region, close to the boundary 
between first- and second-order transitions, the temper- 
ature To is of order BD/{4tt'^), where B is the exponent 
in the expression for the rate. In a typical tunneling 
experiment with a macroscopic lifetime of a metastable 
state, B ~ 477^, so that Tq ~ D. The latter constant can 
be expressed in terms of the anisotropy field and the total 
spin of the particle: D = gpLBHAl{2S). For Ha ^ n 
and S ~ 100 the transition temperature will be of the or- 
der of 10 mK. This is within experimental reach. Notice 
that the smallncss of To in our model comes, in part, from 
the fact that the noncommutation of Sz with the Hamil- 
tonian is small in the lower right corner of the phase 
diagram, where the effect is to be searched for. One 
can expect that in models with transverse anisotropy ob- 
servable first- and second-order transitions will coexist 
at higher temperatures, since the transverse anisotropy, 
rather than the required small transverse field, will drive 
the decay of the metastable state. Such a model requires 
a different approach and will be worked out elsewhere. 
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